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We study the phase diagram of quark matter at hnite temperature (T) and chemical potential (/r) 
in the strong coupling limit of lattice QCD for color SU(3). We derive an analytical expression of 
the effective free energy as a function of T and p., including baryon effects. The finite temperature 
effects are evaluated by integrating over the temporal link variable exactly in the Polyakov gauge 
with anti-periodic boundary condition for fermions. The obtained phase diagram shows the hrst 
and the second order phase transition at low and high temperatures, respectively, and those are 
separated by the tri-critical point in the chiral limit. Baryon has effects to reduce the effective free 
energy and to extend the hadron phase to a larger p direction at low temperatures. 
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I. INTRODUCTION 


Exploring various phases of quark and nuclear matter 
has recently attracted much attention both theoretically 
and experimentally. In the Relativistic Heavy-Ion Col¬ 
lider (RHIC) experiments, it is probable that strongly 
interacting Quark Gluon Plasma is created in heavy-ion 
collisions [j . The phase transition from hadron phase to 
QGP at high temperatures and at zero barwn chemical 
potential is predicted from lattice QCD [1, 0, and vari¬ 
ous experimental signals at RHIC suggest the formation 
of QGP. On the other hand, compressed baryonic mat¬ 
ter is created in heavy-ion collision experiments at lower 
energies, and cold and dense baryonic matter is realized 
in the core of neutron star. For these large baryon den¬ 
sity matter, various interesting matter forms have been 
proposed so far. These states include admixture and su¬ 
perfluidity of hyperons and strange quarks in neutron 
star core [3, the ^P 2 neutron superfluidity [^, pion 
and kaon [7| condensations, color ferromagnetic state [8|, 
color superconductor (GSG) 0, in addition to the for¬ 
mation of baryon rich QGP [lol . ini. 

The lattice QCD Monte-Carlo simulations are possible 
for hot baryon-free nuclear matter, and matter at small 
baryon density [H, [l^l can be studied by, for exaniple, 
the Taylor expansion method around p = 0 [S El, 
analytic continuation method El, canonical ensemble 
method El the improved reweighting method El- 
However, properties of higl^ compressed matter are still 
under debate i El El, El El- This is because the 
fermion determinant, which is used as the weight in the 
Monte-Carlo simulation, becomes complex at finite chem¬ 
ical potential El El- Thus, in order to attack the prob¬ 
lem of compressed baryonic matter, it is necessary to 
invoke some approximations in QCD or to apply some 
effective models . A ]^sible approach is to study 

color SU(2) QCD [12L [^. l22l. [2^, where the fermion 
determinant is still a real number even at finite baryon 
chemical potential, but there are several essential differ¬ 
ences between color SU(2) and SU(3) QCD. For example, 
the color anti-symmetric diquark pair becomes color sin¬ 


glet, whose nature would be very different from those 
discussed in the context of CSC, and this diquark pair 
is nothing but a baryon which is a boson in color SU(2) 
QCD. 


One of the most instructive approximations to inves¬ 
tigate the finite temperature T and chemical potential p 
of QCD is to consider the strong coupling limit of lat¬ 


tice QCD p, l^,M \2' 


In e free en¬ 

ergy at finite T of strong coupling lattice QCD was an¬ 
alytically derived and predicted the second order chiral 
symmetry restoration temperature [H, P . The strong 
coupling limit lattice QCD effective action for finite p 
and zero temperature (T = 0) was also derived with 
the help of lattice chemical potential EB and predicted 
a phase transition near the density of baryonic forma¬ 
tion 13011. These investigations tri gge red many later an¬ 
alytic [Bjp m, [m [m m, H M m, P and semi- 
analytic m investigations of finite T and p of lattice 
QCD in the strong coupling. It is worth to mention that 
the Monte-Carlo numerical results of lattice QCD should 
reproduce the analytic result of the strong coupling, and 
the qualitative nature, and even quantitative nature for 
some physical values such as meson masses, are quite 


close to the reality of finite coupling results [2^ [T 


Based on these past experiences, there have been re¬ 
cently renewed interests of strong coupling lattice QCD 
as an instructive guide to finite T and p QCD. The effec¬ 
tive free energy of strong coupling limit of lattice QCD 
action for SU(2) was analytically derived in [2^. The 
effective free energy of finite p and zero temperature 
(T = 0) for SU(3) was derived by Azcoiti et al. P, who 
developed a method to decompose the coupling term of 
the baryon and three quarks into the coupling terms of 
diquark auxiliary field {4>a) with two quarks and those 
of (fa with a quark and a baryon. The effective free en¬ 
ergy at finite T and p for SU(3) was obtained in [p . 
but the baryon effects are ignored there. Thus there is 
no work which takes account of both finite temperature 
and baryon effects in the strong coupling limit of lattice 
QCD for color SU(3) yet. 
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In this paper, we study the phase diagram of quark 
matter at finite temperature (T) and finite chemical po¬ 
tential (/r) in the strong coupling limit of lattice QCD 
for color SU(3). We derive an analytical expression of 
the effective free energy as a function of T and /r. We 
take account of both the mesonic and baryonic composite 
terms in the 1/d expansion of the lattice QCD action, and 
perform the temporal link variable {Uq) integral exactly 
in the Polyakov gauge with anti-periodic boundary con¬ 
dition for fermions, while we ignore the effects of finite di¬ 
quark condensate. Firstly, our treatment is different from 
the works by Nishida and Bilic et al. who 

extensively studied the phase diagram with the leading 
term in the 1/d expansion containing only the mesonic 
composites. Secondly, our formulation is different from 
the work by Azcoiti et al. [s^, who made the one link 
integral also for Uq which is an approximate treatment at 
finite temperatures. Thirdly, we propose a way to include 
the diquark condensate as a color singlet order parameter 
in Subsec. IIVDI although extensive study is not carried 
out and will be reported elsewhere. 

This paper is organized as follows. In Sec. im we de¬ 
rive an analytical expression of the effective free energy 
in the strong coupling limit of color SU(3) lattice QCD 
with finite temperature and quark chemical potential. In 
sec.nni we study the phase diagram of strong coupling 
limit lattice QCD in the chiral limit. In Sec. CYl we ex¬ 
amine the parameter dependence of the present model 
and compare our results with those in other treatments. 
Also we propose a formulation to include diquark conden¬ 
sates in a mean field ansatz. We stress the importance of 
baryon effects in the phase diagram. We summarize our 
results in Sec. El 


II. EFFECTIVE FREE ENERGY IN THE 
STRONG COUPLING LIMIT OF LATTICE QCD 

In this section, we derive an expression of the effec¬ 
tive free energy in the strong coupling limit lattice QCD 
with Nc = 3 with finite temperature and quark chemi¬ 
cal potential in a mean field ansatz including the baryon 
effects. Chemical potential is introduced in the same 
way as in Ref. [l^. For the finite temperature treat¬ 
ment, we follow the work by Damgaard, Kawamoto and 
Shigemoto [H, [2^ , in which the anti-periodic boundary 
condition for fermions is exactly treated and the integral 
over the temporal link variable Uq is performed exactly 
in the Polyakov gauge. In order to apply this technique, 
we have to obtain the effective action in the bilinear form 
of the^ark field. Such effective actions have been de¬ 
rived [ 2 ^ [ 2 ^ [U [H, m, [H, [ 3 ^ only with the leading 
order mesonic composite term in the 1/d expansion for 
color SU( 3). We utilize the idea proposed by Azcoiti 
et al. [33 to decompose the baryonic composite term. 
Throughout the paper, both of the temporal and spatial 
direction points on the lattice, (3 = 1/T and L, are as¬ 
sumed to be even integers, and the lattice spacing is set 


to be unity. While T = 1//3 takes discrete values, the 
effective free energy is given as a function of T (and /r), 
then we consider T as a continuous valued temperature. 


A. Strong coupling limit and integral over spatial 
link variables 


We start from an expression of lattice QCD action with 
one species of staggered fermion for color SU(Vc). In the 
strong coupling limit {g —> 00 ), we can ignore the pure 
gluonic part of the action, since it is proportional to 1/■ 
As a result, the lattice action contains only those terms 
including fermions, Sp. 


Sf[U,x,x] 


qiUo) 

Jp 



aim) 
O P 


i=i 

\ X! [x{x)e^^Uo{x)x[x + 6 ) 

X 


-Xix + 0)e ^U^{x)xix) 


:'^Vj{x) x(x)Uj(x)x(x+j) 


-x{x + j)U^{x)x{x) , 


( 1 ) 

( 2 ) 

(3) 

(4) 


where we introduce the chemical potential /i in the same 
way to Ref. [3. And ^ (j = 

1, 2, 3) is a Kogut-Susskind factor. The staggered fermion 
X represents the quark field, and the SU(Vc) matrix 17^ 
represents the gauge link variable. 

In the first step, we perform the group integral for 
spatial link variables, Uj{x) {j = 1,2,3). Integral of the 
leading and next-to-leading order terms in the 1/d ex¬ 
pansion leads to the following action, 


= -\{M,VmM) - {B,VbB) . ( 6 ) 


where the inner product of fields are defined as 
{A,VB) = yA{x)V{x,y)B{y). The mesonic and 
baryonic composites and their propagators are de- 
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fined [H, [2^ as 


Dn 


M{x) = x“(a;)x“(a;) , (7) 

B{x) = -^Sab-cX°‘{x)x^{x)---x''{x) , (8) 

B{x) = -^eab-cX‘'{x)---X^{x)x‘'ix) , (9) 

1 

VM{x,y) = > (10) 

i=i 

j=i ^ ^ 

X ('S„+i + (-I)"-'*.,.-;) ^ (11) 

Here we have utilized the S\J{Nc) group integral formu¬ 
lae, 


1 d[U]UabUl, = 

j^r ^ad^bc 1 

» c 

(12) 

1 d[U]UabUcd-- 

‘ ^e/ = ■ 

(13) 


The baryonic composite action {B, VbB) is often ig¬ 
nored with Nc>3, since it is proportional to 
in the 1/d expansion • This scaling can be understood 
as follows. Mesonic and baryonic propagators contains 
the sum over j = 1, 2,... d, and they are considered to 
be proportional to d. In order to keep the mesonic term 
{M, VmM)/2 finite in the large d limit, the mesonic com¬ 
posite should be proportional to . Then the quark 
field, the baryonic composite, and the baryonic composite 
action are proportional to and 

respectively. For the discussion of dense baryonic mat¬ 
ter, however, we expect larger baryon effects. Thus we 
keep this baryonic composite action and proceed. In the 
following discussion, we consider Nc = 3 case. 


B. Auxiliary fields 

The effective action Eq. © contains six fermion terms, 
while we have to obtain the effective action in the bilinear 
or pfafhan form in x to perform the quark field inte¬ 
gral at finite temperature. In order to reduce the power 
in X and x, we introduce several auxiliary fields. 

The highest power term BB containing six quarks can 
be reduced by introducing the auxiliary baryon field b 
through the following identity, 

e(BTBB) ^ detVB J 


Da = ^EabcX^X^" + , DI = '^SabcX^'X^ + ^^X“ ■ 

^. . . ^ 

These are the combinations of diquark and baryon- 

antiquark (antibaryon-quark) pairs, and have the color 
transformation properties of 3 and 3 for Da and D\, 
respectively. The parameter 7 is introduced so as to gen¬ 
erate the coupling terms, Bb + bB. 

D\Da = Bb + bB + Y , (16) 

'y'^ 1 - 

Y = —M^ - -Mbb . (17) 

2 972 ^ ^ 

The decomposition in Ref. corresponds to 7 = 2. 

The product D\Da can be generated by the auxiliary 
field (j)a, and we can replace bB + Bb terms as follows. 

e^B+Sb ^ J ^ (Ig) 

where the expectation value of (pa is the same as that for 
Da, {(pa) = {Da)- 

In terms of the 1/d expansion, the baryonic auxiliary 
field b is proportional to d^^^, provided that the exponent 
in Eq. d is 0{d Thus the second term x& in 

Da is 0 ( 1 ), while the first term SabcX^X^ is proportional 
to and we expect the dominance of the second 

term for large d. This may be the reason why we need 
the baryon-antiquark pair in discussing the diquark pair 
condensate. 

In the next step, we decompose the_ coupling term of 
the baryon and mesonic composite, Mbb, by introducing 
the baryon potential auxiliary field to though the identity, 

^Mbbl9'y^ _ j /2-uj{aM+gu,bb)-a‘^M'^/2 ^ 

where g^j = and (w) = —{aM + ga,bb). 

Note that the local four baryon term becomes zero, 
b{x)b{x)b{x)b{x) = 0, due to the Grassmann variable na¬ 
ture, which is a natural consequence of staggered fermion 
formulation for one-flavor. 

Einally, we introduce the auxiliary field for chiral con¬ 
densate. It is interesting to note that we have additional 
“mass” terms, ( 7 ^ -|- a^)M'^/2 for the mesonic compos¬ 
ite M through the decomposition of baryonic composite 
action by introducing the auxiliary diquark and baryon 
potential fields. These terms are made of four quarks, 
and it is not easy to handle in the quark integral. There¬ 
fore, we include these terms in the hopping term, 

]^{M,VmM) - Y a^)M'^ = ^{M,VmM) , ( 20 ) 

then it becomes possible to bosonize as. 


Next, we decompose the coupling terms of the baryon 
and three quarks by using the technique developed in 
[ 3 ^. We consider the following composite diquark field 


gi(M.yMM) _ J 25j^jg-^(<T,VM<T)-(<T,VMM) ^ ^21) 

VM{x,y) = VM{x,y) - ( 7 ^ +a‘^)Sa:,y . ( 22 ) 
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The expectation value of a is given as (cr) = — (M) . The 
mesonic propagator Vm has negative eigenvalues as well 
as positive ones, and thus it is expected that instabil¬ 
ity is introduced in the gaussian integration. However in 
the mean field ansatz, vacuum expectation value of the 
meson is introduced so that next neighboring {x,x + jl) 
dependence is suppressed, which corresponds to the sup¬ 
pression of (x, y) dependence in Vm{x, y). In this way, we 
circumvent the instability, which is done in the literature 
and we show in the following. 

After these sequential introduction of auxiliary fields, 
we obtain the action of quarks, baryons, diquarks, baryon 
potential, and the chiral condensate as follows. 

Sf = + 4"^ , (23) 

b, (j), ^^cr, w] 

= ib,V^^b) + {4>\(t>) + -f i((T, Vmct) (24) 

S^F [^ 0 , X, X. b, b, (j), (t)\ a, u)] 

= 4“ + (mq, ^ [(X“, i'lb) + {Ha, X“)] 

+ '^Scab [{He, x“x^) + {HH, H)] , (25) 

rriq = Vmo- + auj + mo . (26) 


Fourier transformation for the fermion fields [1^ [2^ , 

e*'="*^V’m(x) , (30) 

1 ^ 

= 7^ E e-^'='">-(x) , (31) 

V ^ m=l 

where ip stands for x or 6 , and the Matsubara frequencies, 
km = 2'iT{m—lj2)jf3, are selected to satisfy anti-periodic 
condition of fermions, '0(/3,x) = —^/’( 0 ,x), to introduce 
temperature effects. 

We ignore the time dependence of bosonic auxiliary 
fields, (p^H (static approximation), and we work in 
the Polyakov gauge, where the link variable Uq is diagonal 
and independent on time, 

t/o(x, r) = diag(e*®^(’')/^, e*® 3 (x)// 3 ) ^ ( 33 ) 

with the condition 9i + 62 + 83 = 0- The quark action is 
found to be represented in the form of pfaffian, 

4 ^ = ^ E (Xm.Xm)G~ 6 ^("^>’^) fe) 

x.m,n ^ 

+ E&- + ^-C'“), (33) 

x,m 


where and Sp'^ are the action of pure auxiliary fields 
and the action containing quarks, respectively. The in¬ 
verse baryonic propagator is modified as 


^ y) =^3^ y) + , 



(27) 


It is noteworthy that the quark action is decomposed 
into that for each spatial point, x. Therefore, it would be 
good enough to assume that bosonic auxiliary fields have 
constant values, i.e. mean field ansatz would be valid. 
This simplifies the term containing Vm’x as follows, 


1, ~ , /3L3 „ 

VmCT = ttaCr , -(cr, VmCt) = 




(28) 

(29) 


C. Quark integral 


In order to perform the quark integral, we would like 
to separate the action into terms, each of which has as 
small number of quark fields as possible. In the quark 
action Eq. (|25]), the time component of the link variable 
C/o connect the quark field of different (imaginary) time, 
and all the quark fields with different times couple. 

This coupling is known to be separated by using the 


Ga{,^(TO,n;6'a) = 
H“(/c) 


^ {kra)bab^mn 'J^cabP^c^mn \ 

^ ^^cahHc^mn ^ {^rriHab^mn J 

(34) 

TUq + i sin(A: -|- 0“/ f3 — iy) , (35) 


= ^Habn. , = ^bmH • (36) 

In the first line of Eq. (1551) . we have used the notation 
5'mn = Sm, 0 -n+i- The second term in Eq. (l33l) can be 
absorbed into the first term by shifting the quark field at 

(C) 

a cost of producing another term Sp , which is bilinear 
in b and b. 



qiC) 

Jp 


^ E Hm,Xm)Gab{'nT'^H 

x..m,n 



+ 54 (,37) 


E {CH-CHGab{ni,n) 

^,m,n 



.(38) 


(C) 

The action Sp appears from the baryon-quark coupling 
generated by the diquark condensate, and it is very diffi¬ 
cult to handle with finite values of p. Another treatment 
to replace this coupling by other terms will be discussed 
in Subsec. lIVDi and we temporarily ignore here. In 
this case, by symmetrizing for m and m' = [3 — m + 1 in 
Eq. (|3^ . We obtain the block diagonal form of Sp, 


/3/2 

4 = E E<^XmWm')gab(fcm) 

X,a,& 771=1 



( 39 ) 
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Sabik) = 

^ab 


^ab 

B{-k)ab 


(40) 


Here, {B{k))ab = 5ahB°‘{k), {A)ab = iecab4’c and we 
have used the relation, = ‘^ir — km = — fcm(niod(27r)). 
Since Xm, Xm' are independent of each other, Grassmann 
integral over XtX leads to a determinant: 


V[x,x]e-^-' 


(41) 


13/2 (3 

G(x) = det [gab(A:™)] = det [gob(fcm)]^^^ • (42) 

m=l m—1 

The G(x) is evaluated by the direct calculation of det: 


G(x) = X{[i\Bi\<tbi\^ + B2\<t32\''+ BM) 

m—1 

+ 7^ E BaB'M{BbB', + B',B,) 

{a^b,c)—cyc. 


+ BiB2B3B[B'^B'3 


(43) 


where Ba = B°‘{krn), B'„ = B°'{—km)- In a similar way 
to that in Ref. [^, we can perform the Matsubara fre¬ 
quency product 


G(x) = n [1 + (44) 

3 

where Zj(x) is the solution of det [gab(fcm)] = 0. The 
explicit derivation is given in the Appendix [XI 


D. Effective free energy at zero diquark condensate 


There are three comments on the phase of G(x) in 
order. First, the square root in the quark determinant 
(Eq. is the Pfafhan root [i^, but it comes from the 
square root in Eq. SI]), where we extend the range of 
the Matsubara product from /3/2 to /3 by using the even 
function nature of det [gab(fcm)] = det [ga6(-^m = fcm')]- 
Since the phase of the square root in Eq. S21) is defined to 
reproduce the product up to /3/2, there is no phase am¬ 
biguity. This is because we can represent the quark ac¬ 
tion S'p ^ in Eq. (I39|) in a usual bilinear fermion action by 
defining a new fermion field as (V'^,, Vil/) = (-xL'. Xm')> 
then it is not necessary to introduce the Pfafhan root. 
Secondly, we may have a phase coming from a constant 
in log G(x) as shown in the Appendix but this con¬ 
stant does not depend on the gluon conhgurations 0a (x:). 
As a result, we have a hxed phase for a given spatial point 
X, and we get well-dehned integral of G(x) over dUo{x.). 
In this way, we expect that we get reasonable continuum 
limit, which we can still consider in the strong coupling 
region. Thirdly, in order to take one flavor fermion con- 
hguration, we have to take one quarter root of Eq. (gSI), 
where the well known phase ambiguity of complex num¬ 
ber appears [4^. Here, we simply consider conhguration 
with one species of staggered fermion, which in general 
regarded as four flavor conhgurations, and do not take 
into account the complexity of the four havor feature of 
staggered fermion. 

By using the SU(3) Haar measure in the Polyakov 
gauge. 


J dUoix) = J , (48) 

A = 6{9i -b 02 + 03 ) J^(I ~ cos(0i — 9j)) , (49) 

i<j 

the integral over Uq can be performed analytically. 


When the diquark condensate is zero, (pa = 0, we know 
the solutions of det [gab{km)] = 0. We can take Ba = 
B°‘{k) = 0 and B'^ = B°‘{—k) = 0, and we get four 
solutions for each a, 

iPza = ±[i{9a-il3/i)±l3Eq{mq)] , (45) 

where Eq{mq) = arcsinbrn^ is one dimensional quark 
excitation energy. We can then explicitly write the quark 
integral results as. 


Nc 


a—1 

1/2 

X (1 + _|_ ^iiSa.-il3lJ.)+l3Eg'^ 

Nc 

= ]/[ 2 (cos(0 — i/3/r)-b coshEg) 

a—1 

= ]/[2(Ga-bG^cos0a-biS'^sin0a) , (46) 

a 

Ca = cosh l3Eq , Cfj, = cosh/3/r , 5^ = sinh /3/i . (47) 


exp ^n/'^G o(x)G(x) , 


= -Tlog 


1 


1 


- ( G;J - -Ga + -Cn, 
Ned 


Ca = cosh ^ , Cn cbi = cosh ^ 


(50) 

(51) 

(52) 


where T = 1//3 is regarded as a temperature. It is inter¬ 
esting to hnd that G(x) can have a complex phase for a 
given gluon conhguration, but after integration over the 
temporal link variables, we have a positive dehnite re¬ 
sults and we do not have the sign problem. This is one of 
the merits in the strong coupling limit in which the link 
integral can be performed in an analytic manner. This 
result is consistent with that for SU(A^c) shown in, for 
example, Ref. [^ . 


mG) fsinh[(iVa + l)A/T] 

- -T log I-3i^h[E/r]-+ 2^"^= 


, (53) 


while CNcbi term does not appear in V{Nc) case [28|, 
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When the diquark condensate is zero, (j)a = 0, we can 
iC) 

ignore Sp , and it becomes possible to perform baryon 
integral, too: 

( 5 a;w) = ^ log Det [1 + gu^ujVB] ■ 

As shown in the Appendix [B1 we can evaluate this de¬ 
terminant by using the Fourier transformation. For large 
spatial lattice size L, by replacing the sum over k by the 
integral, we get the following expression. 




(47rA3/3) Jo 

= -a!,"'/'” ( 


Airk^dk log 




16 


(54) 



a 


where, = 1.0055 , A = 1.01502 x 7t/2, and is 

given as. 


= ilog(l + a:2) 


arctan x — x + — 

O 


• (55) 


Since this baryon determinant term F^g is independent 
from T and /r, it would be more convenient for the later 
discussion to separate the quadratic term in oj as follows, 

Y + +AFjg^(5^w) , (56) 

, 3 (b) /g^Ay 

au; = 1- 5«o ( — I : (57) 



K 


where = 0(0!“^) at small uj values. The second 

term in comes from the expansion of Eq. (1541) . 

After quark, gluon and baryon integral, the total effec¬ 
tive free energy is obtained as, 

FeS = -k (m,) + AF^ (g^uj) , (58) 

where niq = GcO + aw + toq. 


E. Stability and equilibrium condition 

We have introduced two parameters, 7 and a, and two 
auxiliary fields, a and w, in the derivation of the effec¬ 
tive free energy, Eq. (l58l) . Since we have introduced these 
parameters and fields through identities, the final results 
should not depend on these parameters if all the integrals 
are completed. However, we are working in the mean field 
ansatz, then we may have some parameter dependence. 
In principle, we should select the parameters so as to 
keep the mean field ansatz valid; the effective free en¬ 
ergy should be stable against the variation of the fields, 
cr and (u, and the effective free energy (the free-energy 
density) at equilibrium should be stationary against the 
variation of parameters. We further require that the chi¬ 
ral symmetry is restored at very high temperatures. The 


FIG. 1: (Color online) Parameter range which satisfies the 
conditions of stability and high T chiral symmetry restoration. 
The solid dot represents the parameter set, q = 0.2, o? - 1 - 7 ^ = 
1/2 — 0 , which we adopt in the later discussion, and the square 
shows the parameter set without baryon effects. 


stability of the effective free energy against the variation 
of tJ is satisfied when Uo- > 0, and the chiral symmetry 
restoration at very high temperatures is ensured when 
a^j > 0. The region of a and 7 which satisfies both of 
the conditions are shown in the upper panel of Fig. [T] 
The parameter dependence on parameters are discussed 
in Sec. IIV Al 

By using the equilibrium condition, two auxiliary fields 
are related to each other, then we can obtain the effective 
free energy as a function of one order parameter. At 
equilibrium, the effective free energy is stationary with 
respect to a and w, 


dFeS 

da 

dFeS 


aF 


(9) 


= a„a+a„ 


eff 


= a^jW + a 


dmq 


= 0, 


dAF 


(b) 

eff 


dw dniq dw 


= 0 . 


(59) 

(60) 


The effects of AF^g is small when the fields are small. 
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then in this case oj can be represented by the chiral con¬ 
densate a as, 


(9) 


- OJ ~ CT —--- 

a OlTia 


duiQ 


(61) 


With this approximation for 00 , the effective free energy 
is given as a function of a as, 


•T^eff 

rrig 

ba 


(5aCr) , 

b„a + mo , 

v2 


OjQ 


a 

Clfjj 


9<y = 


cyguj 

Uuj 


(62) 

(63) 

(64) 


With this form of the effective free energy, the meaning 
of parameters are a little more clear. The constituent 
quark mass is a linear function of a, then the coefficient 
bcr represents the polarizability of the chiral condensate, 
which is modified by the baryonic composite effects. The 
parameter determines the strength of the coupling of 
the chiral condensate and the baryon, and 
represents the repulsive self-interaction of a coming from 
the baryon integral. 

The parameters 5cr and ga are related to a and 7, and 
they have the region which satisfies the conditions of sta¬ 
bility and high T chiral symmetry restoration as shown in 
the lower panel of Fig.[TJ We notice /S.F^{ga(j) has the 
positive value for any cr, hence, the smaller gn leads to 
smaller Thus, we choose those parameters to give a 
small coupling g^ for a given polarizability bn- The small¬ 
est ga at a fixed ba is obtained in the limit of aa —> +0, 
or 7^-1-^ 1/2 — 0. There is no singular behavior in the 
effective free energy in this limit as a function of the chi¬ 
ral condensate a in Eq. (15^ . For numerical calculations, 
we adopt a = 0.2, which gives almost the smallest ga as 
shown by the filled circle in Fig. [2 as a typical value in 
the later discussion. 



FIG. 2: (Color online) Effective potential as a function of 
a. Upper, middle, and lower panels show the effective free 
energies in the unit of Tc at ^ = 0, T = Ttcp, and 
T = 0, respectively. Dots represent the equilibrium points, 
and dashed lines connect these points. 


A. Zero temperature 


III. PHASE STRUCTURE 

In the previous section, we have demonstrated that 
the effective action in the strong coupling limit lattice 
QCD can be obtained in an almost analytic way with 
1/d expansion and mean field ansatz, when the diquark 
condensate is zero. Especially, we focus our attention to 
the chiral phase transition. 

In this section, we discuss the phase diagram based on 
the effective free energy as a function of the chiral con¬ 
densate in Eq. (ESI) in the chiral limit, toq = 0. Since 
we utilize the linear approximation [a^ojja ~ cr, see 
Eq. (IMTl ). equilibrium value of a slightly differs from 
that of —dFesIdrriQ = —{xx)- However the difference 
of those is small and only by around 1 % as discussed 
in Subsec. IIII Cl For numerical calculations, we adopt a 
parameter choice of a = 0.2 and -|- 7^ ^ 1/2, which 
gives ba — 0.0465 and ga — 2.527. This value of ga is 
almost the smallest value allowed in this model. 


It would be instructive to analyze several limits of the 
effective free energy Eq. (EH). The effective free energy 
from the quark integral F^q depends on the temperature 

and chemical potential. At zero temperature, F^g can 
be reduced to 




-iVc Eg {Eg > g) , 

-W/i {Eg < g) , 


(65) 


where Eg = arcsinhmg is a quark excitation energy. In 
vacuum in the chiral limit, (T,/r, mg) = (0,0,0), 
has a linear term in tr, while other parts of the effec¬ 
tive free energy start with cr^, then we necessarily have 
a finite equilibrium value of cr. On the other hand, for a 
finite chemical potential, F^g becomes independent from 
cr for small cr, and the effective free energy start with the 
quadratic term, ba<j‘^ in the chiral limit. Then we have 
a local minimum at cr = 0 when g is finite even if it is 

































FIG. 3: (Color online) Phase structure in the strong cou¬ 
pling limit lattice QCD with Nc = 3. The solid (red) and up¬ 
per thick dashed (blue) lines show the first and second order 
phase boundary, respectively. The dashed (blue) line shows 
the boundary on which the effective free energy curvature be¬ 
comes zero at a — 0. The dot represents the tri-critical point 
(TCP). Chemical potential and temperature are shown in the 
unit of Tc. 


For larger chemical potentials, fj, > the effec¬ 

tive free energy has a local minimum at cr = 0, as already 
mentioned in the case of T = 0. The above critical chem¬ 
ical potential is shown by the dashed line in Fig.O 

In the present model, the chiral phase transition is sec¬ 
ond order at small chemical potentials in the chiral limit. 
The coefficient C2 is a decreasing function of T for a fixed 
/i for /r/r < 0.588. In addition, the higher order terms 
are positive when the chemical potential is small, 

~ C4a4 + 0(^6) , (71) 

C4 _ 34^^ 20T2-41 + 150/(^3^+ 2) 

bi 28 [AbJ + 12T3(C3^ + 2) 

(72) 

As a result, we do not have any local minimum at finite 
values of a giving smaller effective free energy than ct = 0 . 
Therefore, the above Tc is the actual chiral phase tran¬ 
sition temperature, and then the chiral phase transition 
is second order at zero and small chemical potentials in 
the chiral limit. The behavior of the effective free energy 
at /r = 0 is shown in the upper panel of Fig. 


very small, as shown in the bottom panel of Fig. [2l As 
a result, the first order chiral phase transition occurs at 
the chemical potential which satisfies, 

- (T = 0) = Tesiao;T = 0,/4 = 0) , (66) 

where (Tq stands for the vacuum equilibrium value of a. 


B. Small chemical potentials 

At finite temperatures, we can expand in 

= -T\og 


C: 


3 ^ 


r(C3^ + 2) 


+ a4«) ,(67) 


where ^F^g = 0{a‘^). In the chiral limit, the coefficient 
of (T^ in Fes, 


(rr \ 


( 68 ) 


controls the second order phase transition. At zero chem¬ 
ical potential, this coefficient changes sign at 


Te = = 0) = 


105 „ 


(69) 


For lower temperatures, T < Tc, the coefficient changes 
sign at a chemical potential 


^(2„d)(j.)^|arccosh (^^- 2 ) 


(70) 


C. Phase diagram 

In previous subsections, we have discussed the proper¬ 
ties of the effective free energy at small a in the chiral 
limit. In order to discuss the whole phase diagram, we 
show the results of numerical calculations in the chiral 
limit (too = 0 ) with a parameter value a = 0.2 in this 
subsection. 

Since the chiral phase transition is second order for 
small chemical potentials and is first order for small tem¬ 
peratures, we have to have a tri-critical point (TCP) in 
the phase boundary. At TCP, the finite equilibrium chi¬ 
ral condensate cr, which gives the same effective free en¬ 
ergy as that for cr = 0, approaches to zero. In Fig. [21 we 
show the effective free energy as a function of cr at zero 
chemical potential (upper panel), zero temperature (bot¬ 
tom panel) and at the TCP temperature. We can find 
clear characteristic behavior of the first order phase tran¬ 
sition at zero temperature and the second order phase 
transition at zero chemical potential. At T = Trcp, we 
see a marginal trend. 

In Fig.jSl we show the phase diagram. The dashed line 
shows the critical chemical potential at which the coeffi¬ 
cient of the quadratic term, cr^, becomes zero. Outside 
of this dashed line, we necessarily have a local minimum 
at cr = 0. At low temperatures, we have another local 
minimum at a finite value of cr, giving a lower value of 
the effective potential than that of cr = 0. As a result, 
we have three regions in the (T, fi) plane: The quark- 
gluon plasma (QGP) phase where the chiral symmetry 
is restored, the region of /r < ^i^”'^^(T) where we have 
one local minimum at a finite value of cr, and the region 
^(2nd)(;j^) ^ ^ ^ where we have two local minima. 
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FIG. 4: (Color online) Chiral condensate cr as a function 
of chemical potential and temperature. Thick lines show the 
phase boundary, and the dot indicates the tri-critical point 
(TCP). Chemical potential and temperature are shown in the 
unit of Tc. 


It is interesting to find that, with the current choice of 
the parameter, smoothly decreases as the tem¬ 
perature increases, and it joins with at TCP. In 

the present model with one order parameter cr, the slope 
of /ii^’’*^(T) (i.e. /dT ) in Fig.[3]has to be the same 

as that of in the vicinity of TCP. The first or¬ 

der phase transition condition of equilibrium and balance 
with .7^eff(0) can be solved as 4c2C6 = cl for the effective 
free energy iFeg(cr) = cq -I- C 2cr^ + C4cr^ -I- cecr®. In the 
vicinity of TCP C2,C4 = 0{AT, Afi) are small, then the 
above condition requires C2 = 0{{AT, Afi)^) leading to 
very small C2 which should be on the second order phase 
transition line, provided that cq is finite at around TCP. 
Therefore, negative slope /dT < 0 around TCP is 

a consequence of larger TCP temperature giving a neg¬ 
ative slope of Ttcp > ^ 0.599 Tc. The TCP 

temperature is a solution of a simultaneous equation of 
C 4 = 0 and fj, = /Xc^”'^^(r), 


1 + yi7||^(579^4^^ , 

(73) 

where stands for the first term of C4 in Eq. ((7^ . In 
the present parameterization, the condition Ttcp > 
is satisfied in a wide range, 0.0864 < a < 0.563. On the 
other hand, when we ignore the baryonic action, we get 
cf) = 0, Tc = 5/3 m, then Ttcp - 0.52T^ < T^. 

In Figs.O we show the approximate chiral condensate 
a and the actual chiral condensate (yy) = dTes/dmo 
as a function of the temperature and chemical poten¬ 
tial in the chiral limit. At zero temperature, the effec¬ 
tive free energy is a common function of a in the region 
arcsinh b^a > and then the equilibrium value of a 


Ttcp _ 41 




M/Te 

FIG. 5: (Color online) Approximate chiral condensate cr and 
the chiral condensate —dJ-es/dmo as a function of tempera¬ 
ture (upper panel) and chemical potential (lower panel). Solid 
(red) and dashed (blue) lines show cr and — (xx), respectively. 
Chemical potential and temperature are shown in the unit of 

n. 


stays constant up to p, = As a result, the equilib¬ 
rium free energy is a constant when p, < and de¬ 

creases as Tefi = const. — NcH, leading to a sudden jump 
of the baryon density from = —d!Fes/d{Ncp) = 0 to 
Pb = 1, which is the maximum value on the lattice. This 
behavior is an artifact of the strong coupling limit, and it 
is also found in previous works at zero temperature and 
finite chemical potential [s^, [H, [H, Slj . 


At finite temperatures, a smoothly decreases, and sud¬ 
denly vanishes at pc^***^ when T < Ttcp ■ The chiral con¬ 
densate (Xx) is almost the same as a. This approximate 
relation holds very well for small values of cr, and even 
for large values of a around the vacuum value, the ratio 
changes only by around 1 %. 
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IV. EXTENDED EXAMINATIONS 


The effective free energy derived and examined in the 
previous sections seems to be reasonable, and the cal¬ 
culated results qualitatively agree with those in other 
works. However, there are several unsatisfactory points. 
First, we have to introduce two parameters, 7 and a. 
Several restrictions for these parameters are discussed in 
the previous section, and further discussions is presented 
in Subsec. II V Al Secondly, we hnd quantitative differ¬ 
ences in some thermodynamical variables from those in 
other works. In Subsec. [IV Bl we compare the present ef¬ 
fective free energy and other strong coupling limit models 
proposed so far. Thirdly, while we have shown that the 
diquark effect appears in several aspects of the effective 
free energy indirectly, it is unsatisfactory that we cannot 
treat the diquark condensate directly. In Subsec. IIVDI 
we propose an idea how to include the diquark conden¬ 
sate directly in the effective free energy. 


A. Parameter dependence 


In the previous section, we have shown the relation 
between the scaled variables such as, T/Tc, n/Tc and 
^eslTc- This is because we can remove the major pa¬ 
rameter dependence with these scaled variables at small 
a values. Here we would like to discuss that this scaling 
behavior corresponds to the modihcation of the lattice 
spacing. 

When we explicitly put the spatial and temporal lattice 
spacings (a and at) in the effective free energy, we find 
the following dependence. 




I 

a^at 


\ba{a^cr)^ + AF^g\a^g^a) 

T 

log Gu{ata^baa; atT, atfi) , 


(74) 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


a 

FIG. 6: (Color online) Parameter dependence of /rc(0) in 
vacuum (thick solid line, red), Ttcp (long dashed line, green), 
/ttcp (short dashed line, blue), and a (dotted line, magenta) 
in the upper graph. In the lower graph, thick solid line (red) 
and long dashed line (green) represents the parameter depen¬ 
dence of ba and Qa, respectively. 


where Gu = exp(—F)|g^/T) = / dl7oG(x). The criti¬ 
cal temperature depends on both of at and a as Tc = 
lOaj&cr/Sa^. 

We require that the second order chiral restoration 
temperature should be described independently from pa¬ 
rameter choice, when we choose the temporal and spatial 
lattice spacing appropriately. Actually, the effective free 
energy up to which governs the second order chiral 
restoration at small /r, is found to be independent from 
the parameter choice for a given Tc, 


T 

a- 


— -h /3 cAT^^^ 
20 




-L] r — 

Tc ^ ViO/3c’/9cTc’/3cTc 


Tc 


3^2 T , ^ ■ 

-log Gtj 

20 Tc 


+ 0{a^), (75) 


where a a = a'^ajat stands for the chiral condensate mea¬ 
sured in the unit of atja^ and /3c = I/otTc denotes the 
temporal lattice size at the critical temperature T = T^. 
Since starts from cr"'’ and Gjj is a function of 

arcsinh {baa)/T and fi/T, the scaled effective free energy 
Teff/Tc is a function of scaled variables Ua, T/Tc and 
/i/Tc when we ignore 0 {a^). 

The remaining parameter dependence may come from 
the mean field ansatz. Thus the parameter should be 
chosen in the range where the mean field ansatz is valid; 
i.e. the dependence of obtained quantities is small. In 
Fig. [51 we show the parameter dependence of a in vac¬ 
uum, Trcp, /iTCP and /rc(0), which suffer from higher 
order contributions of a. Most of these quantities have 
extrema at around a = 0.2 (a ~ 0.188 for Ttcpj Mtcp 
and mUc(O), and a ~ 0.193 for a in vacuum). In ad¬ 
dition, we find that the parameter dependence is not 
strong in the parameter range, 0.1 < a < 0.6. It is 
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FIG. 7: (Color online) Parameter dependence of the effective 
free energy in vacuum in the chiral limit, (T, /r, mo) = (0,0, 0). 
Long dashed (green), thick solid (red), dash-dotted (light 
blue), and dotted (magenta) lines show the effective free en¬ 
ergy in the present work, with parameters a = 0.1, 0.2, 

0.5 and 0.6, respectively. Short dashed line (blue) shows the 

('I'') 

effective free energy without baryon effects, 



^i/T, 

FIG. 8: (Color online) Parameter dependence of the phase 
boundary. Long dashed (green), thick solid (red), dash- 
dotted (light-blue), and dotted (magenta) lines show the 
phase boundary in jFefj with a =0.1, 0.2, 0.5 and 0.6, re¬ 
spectively. Short dashed line (blue) shows the boundary in 
^, without baryon effects. Thin solid line (black) indicates 
C 2 = 0, where C 2 is the quadratic coefficient of a. 


worth to mention here that is small enough at around 
a ~ 0.2 and it satisfies the even integer condition for 
Pc = ^jaTc = 3/105o. > 2 in a symmetric lattice. 

In Figs. [7] and [SI we show the parameter dependence of 
the effective free energy in vacuum and the phase bound¬ 
ary. In these figures, we also plot the results with the 
effective free energy , 


_ 

■^eff — 


-b^pa^+F^^Hm,)=Tc 



^logGc/ , 
(76) 


nig = b'^pa + mo, bf'> = Tc = 


m. 


(T) 


2 Nc 


= 3 , (77) 


which is obtained by ignoring the baryon effects and in- 
tegratin g ov er Uq exactly in a similar way to that in 
Ref. [^, [3^ [3^ . It is clearly seen that large energy gain 
is obtained with a ~ 0.2, and the phase boundary ex¬ 
tends to the larger fi direction. When we ignore baryon 
effects, the effective free energy and phase boundary with 
roughly corresponds to those with a = 0.6 in the 
present model with baryons. 


B. Comparison with other treatments 

While we have treated the time-like link variable Uq ex¬ 
actly in the previous section, the anti-periodic boundary 
condition may not be very important when the tempo¬ 
rary lattice size, /3 = 1/T, is very large. In this case, it is 
possible to perform the one link integral also for 3^°'^ as 


other spatial action, = 1,2,3). After introducing 

auxiliary fields, 6, b and cr, we obtain the action, 

J ^ J , (78) 

SF\x,Zb,b,a] = ^'^<7{x)VMo{x,y)cr{y) 

x,y 

+ Y1 [Hx)Vgi^{x,y)b{y) + a{x)VMo{x,y)M{y) 

x,y 

-'^[Hx)B{x) + B{x)b{x)] , (79) 

X 

1 ^ 

^Moix, y) = ^ ( (by^x+C' + by^x-v) , (80) 

■= i'=0 

VBy{x,v) = Vb-\ {e^^5^,x+b - (81) 

In the above action, quark fields completely decouples 
in each space-time point, then it is possible to perform 
quark integral. For example, when we ignore the baryon 
effects and carry out the quark integral, we obtain the 
following effective free energy [^ . 

- Nc log(5i0V + mo) , (82) 

where = {d+ l)/2Nc and m, = b^^a + mg. The di¬ 
verging behavior at cr = 0 in the chiral limit is suppressed 
when we include the baryon effects in a similar way to 
that in [^ . 

T, yi) . (83) 
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FIG. 9: (Color online) Comparison of effective free energies 
in different treatments. Thin solid (red), dotted (magenta), 
long-dashed (green), and short-dashed (blue) lines show the 
effective free energies respec¬ 
tively. For (long-dashed, green), we show the results 

with two parameters; a = = 0.2 (thick) and 7 = 2 

(thin). The thick solid line (red) indicates the effective free 
energy in a finite temperature treatment with baryon effects 
, present work). We show the results in vacuum in the 
chiral limit, (T, /r, mo) = (0, 0, 0). 

The expression of the baryon integral, is shown in 

the Appendix[Bl It is also possible to obtain the effective 
action with diquark field as, 


reported to lie at f = 0 [s^, and the results with two 
parameters are compared; the same parameter set as in 
the present work, 7 = — oP- with a = 0.2, and the 

value originally adopted in Ref. [s^, 7 = 2. From a com¬ 
parison of iFgg and the main role of baryons is 

found to reduce the effective free energy at small <7 val¬ 
ues, in addition to suppressing the diverging behavior at 
cr = 0. On the other hand, in zero temperature treat¬ 
ments with baryons, and , we find a bump at 

a small A, which separates two local minima. This bump 
comes from the slow startup of the baryon contribution 
proportional to cr® at small a in , and from a cancel¬ 
lation between — log 0 and (m) in F^^^^. In a finite 

temperature treatment, F^g , reduction effect is smooth, 
and we find only one local minimum. The effective free 
energy in the present work, = FeS in Eq- is 

smaller than those in other treatments except for F^^^^ 
with 7 = 2, with which the effective free energy becomes 
unstable in a finite temperature treatment. The large en¬ 
ergy gain in may partly come from the scaling of 

Feslbi^, since the 60- in this work is the smallest among 
the models compared here. In order to compare the ab¬ 
solute values of the effective free energy more seriously, 
it would be necessary to fix the lattice spacing and thus 
the energy scale and to include the effects of the higher 
order contribution in the 1 /d expansion. 


C. Expected evolution of the phase diagram 


.7-(0bv) 


0 = 

m = 


+ ^v"] , (85) 

( 86 ) 


1 


Rvi^ ' 9' 

Amq {3y'^/Ry - niq^) 
0 


where Ry = 1 — u^/3. This effective free energy is es¬ 
sentially the same as that in Ref. [s^, while we use a 
different notation and introduce a parameter 7 as in the 
previous section. 

These effective free energies have similar asymptotic 
behavior for large cr in vacuum. Knowing the asymptotic 
form, F^^\in) — log2TO at m ^ 00, we find all the 

potential terms in and , have the form of 

—Nc log a + const, in the large a limit. The effective free 
energy at finite T, F/f^ in Eq. (1751) . also has the potential 
of the above form, since F^g (a) — log2cr(cr ^ 00). 

In Fig. 1 we compare the effective free energies as 
a function of the chiral condensate a in vacuum in the 
chiral limit. We show the scaled effective free energies 
F^g/ba'^ instead of F^g/Ty, since the chiral restoration 
does not emerge with zero temperature effective free en¬ 
ergies, , F ^^^, • In , the diquark conden¬ 

sate V is set to be zero, as the global minimum is already 


One of the common problems in the strong coupling 
limit is the too small critical chemical potential fic{T = 0) 
relative to Ty = Ty{fi = 0). In Table HI we compare 
the ratio of the critical baryon chemical potential at zero 
temperature with respect to the critical temperature at 
zero chemical potential, = 3/ic(T = 0)/rc. All the 
models based on the strong coupling limit give much 
smaller values for this ratio, < 1-5, than the em¬ 
pirical value, = (1 — 2)GeV/170MeV ~ (6 — 12). 
In the Monte-Carlo simulations at finite quark chemical 
potentials with finite 1/5^, it is not yet possible to ob¬ 
tain fiy{T = 0), but larger values are suggested. 
For example, several Monte-Carlo methods are in agree¬ 
ment with each other for small quark chemical poten¬ 
tials Ai/r < 1 [H, [11 [H Hi, HI , and the critical tem¬ 
perature for these chemical potentials are large enough, 
Ty{^)/Ty{fi = 0) > 0.9, implying that 3. 

Thus for a quantitative discussion, the strong coupling 
limit in the chiral limit with one species of staggered 
fermion is not enough, and it is necessary to take care 
of finite quark mass toq, multi-staggered fermions, finite 
l/g^, other order parameters than the chiral conden¬ 
sate, and/or other mechanisms towards the real world 
in order to explain large Hy{T = 0) relative to Ty, as 
illustrated in Fig. [TOl With finite quark mass toq, the 
effective free energy Fyft always has a minimum at fi- 










13 


TABLE I: Ratio of the critical chemical potential at zero tem¬ 
perature fJ,c{T = 0) and the critical temperature Tc = Tc{i-i = 
0) in strong coupling models. In we have assumed 

Tc — 5/3 to obtain the ratio. For the results in Ref. [3^ . 
values of critical lattice anisotropy a/at and that multiplied 
by are taken for the number of staggered fermions, / = 1, 2 
and 3. 


Model 

Tc 

^^c{T = G) 

3fic{0)/Tc 

^ [30] 

^ m 

^(Tb) 

5/3 
10b,,/3 

0.66 

0.33 Tc 

0.45 Tc 

1.19 {Tc = 5/3) 
0.99 

1.34 (a = 0.2) 

Ref. [3^ 

2.57 

0.57 

0.67 (/ = !) 


2.19 

0.57 

0.78 (/ = 2) 


2.07 

0.57 

0.83 (/ = 3) 

Empirical 

170 MeV (l-2)/3 GeV 

^ (6-12) 


(see Fig. [S]) have to be the same as that of at 

the vicinity of TCP in the chiral limit with one order 
parameter. For this point, there is a debate between 
two Monte-Carlo simulations, one of which suggests the 
smooth connection of cross over boundary and the first 
order boundary [3, [TsI . IT^ , and the other suggests a fi¬ 
nite difference in slope [l7j| . Provided that the nature of 
TCP remains in CEP even with finite quark mass, our 
discussion in Subsec. UlI Cl supports the former if there is 
only one order parameter. However, both of the above 
two results can be consistent if there are other order pa¬ 
rameters than the chiral condensate. Smooth connection 
is expected for chiral transition in methods based on the 
analyticity while we may see other transi¬ 

tion in a direct Monte-Carlo method at finite chemical 

potentials [i3- 



FIG. 10: (Color online) Expected phase diagram evolution 
from the strong coupling limit in the chiral limit with one 
species of staggered fermion towards the real world. 


D. Color angle average in diquark condensate 

In deriving the effective free energy in Eq. HMD we have 
assumed that the diquark condensate takes zero values, 
(pa = 0- If we ignore the diquark-gluon-baryon coupling, 
Sp in Eq. (l38l) . it is possible to obtain the solution of 
detgab(fcm) = 0 and the integral over Uq numerically. 

(C) 

With sp , however, since Uq depends on x and baryon 
fields are spatially connected through Vb, we have to 
carry out the integral of baryon determinants over 2 L^ 
dimensional variables in f7o(x). In this subsection, we 
would like to show an idea to solve this problem in a case 
with one species of staggered fermion. 

Since the diquark field (pa is not color singlet, its aver¬ 
age over the color space should be zero. Thus we cannot 
treat it as an order parameter. One of the way to rem¬ 
edy this problem is to carry out the integral of the ” color 
angle” variables in (pa, then only the color singlet com¬ 
bination = (pl^cpa remains. It is possible to carry out 
this color angle average in a straightforward way. 


nite cr, then the second order boundary becomes cross 
over and the tri-critical point (TCP) becomes the criti¬ 
cal end point (CEP). In addition, since the finite quark 
mass rriQ increases the baryon mass which is closely re¬ 
lated to Tc{T = 0) d^, finite mg is believed to increase 
Hc{T = 0), as shown for example in Ref. [s^. With 
multi-staggered fermions, is suppressed as discussed 
in Ref. |34l |. It would be natural to expect that de¬ 
creases as I/g^ grows, because hadrons and glueballs are 
more bound at larger couplings and thus hadronic phase 
would be the most stable in the strong coupling limit. 
We further expect that the finite coupling effects appears 
most strongly at /i = 0, where the role of gluons relative 
to quarks is the largest. Actually in Ref. [33, it is shown 
that T'c(m) decreases as I/g^ increases, and this reduction 
is more rapid at /r = 0 than at finite fj,. 

With other order parameters than the chiral conden¬ 
sate, the phase diagram will have a richer structure. In 
Subsec. IIII Cl we have discussed that the slope of 



Here we explicitly show the sum or product over color 
indices, and (•••)« means the color angle average. When 
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we integrate over the phase variable for each (j)a, we only 
have those terms having the same power of (pa and (p\ 
as shown in the second line. The power of the diquark 
composite Da and is limited to four as shown in the 
third line by the Grassmann nature. For example, the 
power four terms such as 

= {dId^){dId,) = (88) 

already contain all the Grassmann variables, and product 
with other Da or vanishes. By using Eq. ((551) . we find 
that the fourth order terms in Da and can be arranged 
in the form of (X)In the second order 
terms, we use the symmetry for each color index, for 
example, {(p\pi)v = u^/3. 

Unfortunately we again have the term containing the 
coupling of three-quark and baryon, bB Bb, from 
DlDa = Y + bB + Bb (see Eq. ([I2l)). 


exp(6B -P Bb) 


= / V[v\e ” ‘ ( e 




= J 'j:>)yj(,-^^-^+'"^(bB+Bb+Y)/3+v*M^bb/162 

^ ^v^{bB+Bb)/3-v^-R^Y+v*M^bb/162 


Ry = I — v'^/3 . 


(90) 


In the third line, we have assumed that the integral in 
r.h.s. can be approximated by the representative value of 
V. This approximation would be valid when the diquark 
condensate is strong. In this mean field ansatz, then we 
can solve this self-consistent relation as follows, 

^R^(bB+Bb) ^ ^-v^-RyY+v^MHb/162 /gj^N 


where we have ignored the constant shift in the exponent. 
As a result, we obtain the following relation. 


exp{bB -I- Bb) 


exp 


v'^M^bb 

Ry 162Ry 


.(92) 


Coupling terms of M^bb can be bosonized by introducing 
n bosons, whose expectation values are related to the 
product of qq pair and bb pair M^bb, {ujk} = —o.k{M) -|- 
Pk{M^bb), in a similar way to that in Sec. dll After 
introducing three auxiliary fields, a;2, wi, wq, it is possible 
to carry out the Grassmann variable integral, b and y, 
and we obtain the effective free energy. 


r{Tbv) 


Fx 

Qj(j 

rriq 


9uj 


Fx{cr,v,UJ^) + F^g\ga,uj) + F^^ {rriq) ,(93) 


lx 9 9 9 9 >. V 

-(a^cT uj ^ 2 ) ~s~ 

Z rCy 

1 2 2 2 2 

--7 -a -a^-a2 , 


OcrCT + aul F aiUJi -I- a2a;2 -I- mo 
1 


1 + 


2 4 

'y V UJ1OJ2 


9q! 7^ [ 18aiQ;2i?t, 


(94) 

(95) 

(96) 

(97) 


Here we have replaced ujq = to and ao = a- 

With zero diquark condensate, the effective free en¬ 
ergy in Eq. (IMl) . has a similar structure to FeS 

in Eq. ((62l) . Specifically, when we take the linear ap¬ 
proximation in the same way to that in Eq. dSD, we find 
that the chiral condensate polarizability and the coupling 
constant are the same as before, = ba-,g^^^^ = go- 

defined in Eq. (IMll , and we obtain the same effective free 
energy as before defined in Eq. (EH), 

= 0 ) 

^ + mo) + a) 

= F,b . (98) 


This equivalence may serve a cross check of the effective 
free energy derived in Sec. HD 

On the other hand, there is no potential effects propor¬ 
tional to while we have quadratic term in Fx , then we 
will not have any second order phase transition to the di¬ 
quark condensed state. This comes from the cancellation 
in Eq. (1^ in the case of one staggered fermion. De¬ 
tailed analysis of the effective free energy E g. (l93l) and 


its extension with multi-staggered fermions [2£ 
be reported elsewhere. 


will 


V. SUMMARY 

In this work, we have studied the phase diagram of 
QGD for color SU(3) at finite temperature (T) and finite 
chemical potential (/r) by using an effective free energy 
derived in the strong coupling limit including baryon ef¬ 
fects. We have adopted the effective action up to the 
next-to-leading order of the 1/d expansion (0(1) and 
0{l/'/d)), and by using the mean field ansatz, an an¬ 
alytical expression of the effective free energy is derived. 
The baryonic composite term in the effective action is 
decomposed into the terms consisting diquark conden¬ 
sates, baryons, and quarks [s^. By introducing auxiliary 
fields of the baryon, diquark, baryon potential, and chi¬ 
ral condensate, we have obtained the effective action in 
the bilinear form of fermions. Then the Grassmann inte¬ 
gral of quarks and the sum over the Matsubara frequency 
can be carried out exactly, provided that the solution of 
det gab(^m) = 0 is obtained. At zero diquark condensate, 
we can further perform the integral over the temporal link 
variables and baryon fields analytically. 

This is the first trial which introduces baryon and finite 
temperature effects simultaneously in the strong coupling 
limit of lattice QGD for color SU(3). It is important 
to note that baryon has effects to reduce the effective 
free energy FeS as shown in Fig. [7] and to extend the 
hadron phase to a larger fi direction at low temperatures 
as shown in Fig. [51 when Fes, b'- S'lid T are measured 
relative to Tc- We may expect that this feature remains 
in the realistic parameter region of finite 1/g^- It would 
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then be interesting to compare the phase boundary be¬ 
havior between SU(3) and U(3) to examine the baryon 
effects in this parameter region. 

The obtained phase diagram have the second order 
phase boundary at small chemical potentials, and the 
first order phase boundary at small temperatures sepa¬ 
rated by a tri-critical point. This feature is the same 
as that in previous works, but the ratio of the critical 
baryon chemical potential at zero temperature with re¬ 
spect to the critical temperature at zero chemical po¬ 
tential, R^t = = 0 )/Tc{n = 0), is found to be 

much smaller than the empirical value or that suggested 
in Monte-Carlo simulations. Small R^t is a common fea¬ 
ture in models based on the strong coupling limit, and 
it would be necessary to extend in the direction of the 
reality axes in Fig.jTUlfor a quantitative discussions. On 
the other hand, we expect that Monte-Carlo simulations 
should reproduce the strong coupling results of the phase 
boundary including the small value of R^t at a large 
value of g. 

Finally, we have proposed a method, color angle aver¬ 
age in colored auxiliary fields, which enables us to extract 
a color singlet order parameter and to include the diquark 
condensate explicitly in the effective free energy. 

One of the problems which we have found in this work 
is the parameter dependence of the effective free energy 
J-eS- During the bosonization, we have introduced two 
parameters, a and 7. Since these are introduced through 
identities, the results should not strongly depend on the 
parameter choice and in fact we have shown that we can 
absorb a major parameter dependence of at small a 
values, which determines the second order chiral transi¬ 
tion, in the choice of the lattice spacing. For the remain¬ 
ing parameter dependence, we have required that the 
scaled effective free energy TesjTc in vacuum becomes 
as small as possible, and we have adopted a = 0.2 and 
-I- 7 ^ = 1/2 — 0. This choice of parameters results in 
the temporal lattice spacing of Pc = ^/aTc{fi = 0) ~ 6.45 
in a symmetric lattice. However, we have learned that 
l/aTc{g, = 0) is not large and less than two in the strong 
coupling Monte-Carlo simulations with one species of 
staggered fermion (without quarter root of the quark de¬ 
terminant as in the present work) (^ . It means that 
the parameter region, 0.5 < a < 0.6, is preferred, rather 
than a ~ 0.2. Considering these situations, we have to 
agree that the baryonic effect on the phase diagram del¬ 
icately depends on the choice of the parameter a, and it 
is desired to find a general procedure to determine them. 

There are several future issues to be studied further: 
The first one is an extensive analysis of the effective free 
energy at finite quark masses and/or with diquark con¬ 
densates. Secondly, interesting and promising direction 
is to consider multi species of staggered fermions, since 
color superconductor is expected to emerge at high den¬ 
sities when multi quark flavors are introduced [^. 
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APPENDIX A: SUM OVER THE MATSUBARA 
FREQUENCIES 


The quark determinant G(x) appeared in Eq. (I43|l is 
an even function of km, then it may be expressed as a 
polynomial of cos km ■ 

/3/2 0 

G(x) = det [gabikm)] = det [gab{km)]^^‘^ 

m—1 m—1 

0 6 

m—l j=l 

where r = 1. The derivative by r reads 


dlogG(x) 1 


-V. 


dr 


2 ^ cos km — rYi 

m,j 


= —y 

2 Q 


dz —Yj 


—ip 


E 


27rz cos z — rYj 1 -|- 
—ip 


-sinz!: l + 


ip Yj 


1 


2il 2 ^ sinz)' 1 -I- exp(z/3z)') 


-ip dz/ 

— E^ 


20 ^ dr 1 + exp(z/3z/ 
d 1 


dr 20 


^log(l-Hexp(-i/3z/)) . (Al) 




In the contour integral, we have contributions from z = 
km poles as well as z = z/, whose sum becomes zero. 
Here zj is the solution of cosz = rYj, and O stands for 
the degeneracy for z/(- |-2n7rz). In the second line from 
the bottom in Eq. I|A1(1 . we have used the relation, 


d cos z/ dz/ 


(A2) 


Now we obtain logG up to a factor. 

logG(x) = JEE log (1 -I- exp(—z/3zj)) -I- const. . 

(A3) 
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The sum over Zj is understood as we ignore the degen¬ 
eracy 2 mTi, but we still have two solutions in a pairwise 
way, ±Zj, since Zj is the solution of cos z = Yj. We choose 
one of them as a principle value. 

logG(x) = i ^log [(1-h + const. 

3 

= i ^ log(l -I- cos Pzj) + const. . (A4) 

3 


We ignore the constant terms in logG(x), and we get 
G(x) as follows. 


G(x) 


1/2 


n(l + cos/3zj(x)) 

3 


(AS) 


APPENDIX B: BARYON INTEGRAL 


In this appendix, we show how to obtain the baryon 
determinant Det(l -I- ujVb)- 

First, we make a Fourier transformation of baryon 
field, 


&m(x) 


k=— (fci,fc2,fc3) .(Bl) 


The staggered factor r]j{x) in Vb connects four different 
momenta, 

, k(^) = (fci-k TT, fc2, fc3) , 

k(^) = (fci-I-TT, A :2 + tt, fc 3 ) , = (fci,A :2 +7r, A: 3 ) , 

(B2) 

and two different frequencies, m and m+/3/2. Asa result, 
Vb is found to be block diagonal. 


E 




,k,k' 


L/2 L 0/2 

= E E E 

ki,k2 = l ks m=l 





where represents the baryon field with four different 
momenta. 


It is interesting to find that the square of S becomes a 
c-number. 


1 ^ 

S-S = —— s^l, (s^ = sin^ kj) . 


(B6) 




Now we can evaluate the baryon integral, 

exp(—) = DetVs J I?[6,6]exp — /b,Vg^bj 
= Det [1 -I- ujVb] 

L/2 L/2 L 0/2 

^ n n n 

ki — l k 2 — l ks — l m—1 
L/2 L/2 L 0/2 

= n n n n 

ki — 1 k 2 — l ks — l m—1 

= n(i+^v/i6)^/" 

k 


1 wS 
wS 1 , 


(B7) 


For very large spatial lattice size L, we can replace the 
sum by the integral, 




2L3 

k 

1 

= / dh 


1-k 
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I.Tr/2 


271^ 
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J —7r/2 
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/ dks log 
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16 
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Airk^dk log 


/o 


21,2 3 


1 + 


LU'^k 


16 




(B8) 


In the fourth line, we have made an approximation to 
replace the average in a box to that in a sphere. With 
this approximation, the effective free energy from baryon 
integral can be represented by a function, 


^ J k‘^dklog{l + fc^) 


= ilog(l + a;2)-^ 


..31 


arctanx — x 


■ (B9) 


bm ; ^mk(2) 5 5 ^mkC^)) 5 

Bm “ '^m-\- 0/2 ■ (B4) 

The matrix S represents how these different momentum 
states are connected through r\j, 


/ sin ki 

sin k 2 

sin ks 

0 \ 

sin ^2 

— sin ki 

0 

sin ^3 

sin k^ 

0 

— sin ki 

— sin ^2 

V 0 

sin ks 

— sin /c 2 

sin ki / 


(B5) 


From numerical studies, following normalization factor 
and the cut off, 

a/0 = 1.0055 , A = I X 1.01502 , (BIO) 

are found to give a good global fit of . 

The baryon determinant in Eqs. (1551) and (1551) 
can be also evaluated by using a similar technique shown 
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here. We show only the results here. 
= -^ logDet (m/4 + Vb/x) 


2L3 


log [m^ + sin^(fco — 3i/r) + s^] 


fco,k 


= — —g log cosh(/3arcsinh y/s^ _|_ ^ 2 '^ _|_ 


(Bll) 


where, numerical calculations 


in Subsec. IIVBI we have adopted the following approx¬ 
imation assuming that the spatial lattice size L is large 
enough and that the average in a cubic box can be well 
approximated by the average in a sphere, 


/ ^klog[C6(k,m)+C3j(pi2) 


C{, (k, m) = cosh /3arcsinh \/k'^ + m'^ , (B13) 


where the simple ansatz = 1 and A = TTj2 are used 
in the numerical calculations shown in Subsec IIVBI 
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